I 


AD-A015  201 


OPTIMUM  DESIGN  OF  STIFFENED  SHEAH  WEBS  WITH 
SUPPLEMENTARY  SKIN  STABILIZATION 

D.  M.  Richards 

Cranfield  Institute  of  Technology  (England) 


Prepared  for: 

European  Office  of  Aerospace  Research  and 
Development 


25  August  1975 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


ADA015201 


* 0*  lU:  - TR-  7 5 **  f 9 

- 281094 

Contract/Grant  amber:  AFOSR  74-2693 

OPTIMUM  DESIGN  OF  STIFFENED  SHEAR  WEBS 

WITH  SUPPLEMENTARY  SKIN  STABILISATION. 
* 


D.M.  Richards 

Aircraft  Design  Division 

College  of  Aeronautics 

Cranfield  Institute  of  Technology 

Cranfield,  Bedford  MK43  OAL 

England. 

25th  August  1975 


♦ 


Final  Repoit,  1st  Feb.  1974-31st  Jan.  1975. 
(Scientific  Report  No.l) 


Approved  for  public  release; 

Distribution  unlimited 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

* US  C*p»*««**i  ol  Co*»*»u* 

Prepared  for  v*  »•»» 

CRANFIELD  INSTITUTE  OF  TECHNOLOGY 

AllO 

USAF  EUROPEAN  OFFICE  OF  AEROSPACE 
RESEARCH  AND  DEVELOPMENT 
London . England . 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  Report  Number 


2.  Govt  Accession  No. I 3.  Recipient  * Catalog  Nunber 


4.  Title  (and  Subtitle) 

OPTIMUM  DESIGN  OF  STIFFENED  SHEAR  WEBS 
WITH  SUPPLEMENTARY  SKIN  STABILISATION 


5.  Type  of  Report  & Period  Covered 

Final  Scientific  Report 
1 Feb  74  - 31  Jan  75 


6.  Per forming  Org.  Report  Nunber 


7.  Author (s) 


9.  Performing  Organization  Nome  and  Addreaa 

Cranfleld  Institute  of  Technology 
Department  of  Aircraft  Design 
Cranfleld,  Bedford,  U.K. 


11.  Controlling  Office  Naiae  and  Addreaa 

EOARD  (BOX  14) 

FPO  N.Y.  09510 


14.  Monitoring  Agency  Name  and  Address 

EOARD  (BOX  14) 

FPO  N.Y.  09510 


8.  Contract  or  Grant  Number 
AFOSR  74-2693 


10.  Program  Element,  Project,  Task 
Area  & Work  Unit  Nunber s 

61102F 

967-99-35 


12.  Report  Date 
25  August  1975 


13.  Nunber  of  Pages 
26 


15. 


16.  fx  17.  Distribution  Statement 

Approved  for  public  release;  distribution  unlimited. 


18.  Supplementary  Notes 


19.  Key  Words 


Optimal  design,  shear  webs,  buckling  strength,  plate  buckling 


20.  Abstract 

This  report  documents  the  result  of  an  investigation  of  the  optimum  design  of 
shear  web  structures  whose  function  is  to  resist  buckling.  The  structure 
considered  was  comprised  of  a simply  supported  thin  metal  web  of  constant  depth, 
but  infinite  length,  stabilized  at  uniform  intervals  by  vertical  stiffeners. 

The  analysis  was  based  on  that  of  Stein  and  Fralich,  and  Cook  and  Rockey,  and 
takes  proper  account  of  individual  stiffener  action  in  the  buckling  process  and 
utilizes  a realistic  buckling  mode  between  stiffeners. 

Results  Indicated  that  structural  weight  decreases  continuously  as  stiffener 
spacing  was  reduced  until  an  absolute  minimum  was  achieved  at  a spacing  well 
below  practical  limits.  Computer  routines  based  on  the  analysis  have  been 
written. 
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SUMMARY 

The  optimum  design  of  stiffened  shear  webs  is  considered. 

A design  procedure  is  developed,  based  on  the  analysis 
of  Stein  and  Fralich* , supplemented  by  an  orthotropic  plate 
solution  due  to  Timoshenko11. 

The  effect  of  adding  a layer  of  low  density  material  to 
the  skin  between  stiffeners  is  included. 

The  results  show  that  the  lightest  designs  require  very 
small  stiffener  spacing. 

Substantial  improvements  in  efficiency  are  indicated 
when  supplementary  skin  stabilisation  is  utilised. 
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1 . Introduction 

The  problem  of  designing  shear  webs  has  been  studied 

0 

by  many  authors  including  notably  Kuhn) * 2 Symonds)  Shanley* 
and  Rothwell*  and  more  recently,  Parsons  and  Beard*,  and 
Laakso  and  Strayer7.  However  it  has  rarely  been  possible 
to  generalise  the  results  obtained  sufficiently  to  allow 
direct  application  by  the  structural  designer. 

This  report  is  concerned  with  the  optimum  design  of 
shear  webs  of  the  type  shown  in  figure  1. 

The  structure,  which  is  required  to  resist  buckling 
under  the  design  loading,  comprises  a simply  supported  thin 
metal  web  of  constant  depth  hut  infinite  length,  stabilised 
at  uniform  intervals  by  vertical  stiffeners.  The  web’ between 
the  stiffeners  may  be  further  stabilised  by  the  addition  of 
a layef*  of  rigid  organic  foam.  The  structure  is  assumed  to 
be  linearly  elastic  up  to  the  design  load. 

The  analysis  of  buckling  is  based  on  that  of  Stein  and 
Fralich*  and  Cook  and  Rockey’  which  is  incorporated  in  the 
Engineering  Sciences  Data  Sheets10  widely  used  by  structural 
analysts  and  designers.  This  analysis  takes  proper  account 
of  individual  stiffener  action  in  the  buckling  process,  and 
utilises  a realistic  buckling  mode  between  stiffeners. 

Comparison  has  been  made  with  designs  based  on  Timoshenko's11 
orthotropic  plate  buckling  theory  which  assumes  the  effect  of 
stiffener  stiffness  to  be  smeared  uniformly  over  the  web. 

This  method  is  fo«.  id  to,  be  most  appropriate  as  stiffener  spacing 

tends  to  zero,  whereas  buckling  strength  will  be  over  estimated 
progressively  as  stiffener  spacing  increases. 

The  work  reported  here  was  supported  in  part  by  the 
United  States  Air  Force  European  Office  of  Aerospace  Research 
and  Development  under  grant  number  AFOSR-4-2693. 
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2.  Analysis 

2 . 1 Buckling  analysis  with  finite  stiffener  spacing 

An  infinitely  long  shear  web  of  depth  a*,  thickness  t* 
is  to  be  designed  to  resist  an  applied  shear  flow  q*  without 
buckling. 

The  web  is  subdivided  into  panels  by  a series  of 
vertical  stiffeners  spaced  a distance  c*  apart  along  the 
web. 

The  buckling  of  the  web  may  be  characterised  in  terms 
of  Individual  panel  dimensions  by  the  following  relationship 

ta? 

q*  ..  - KE*  ~ ■ /1X 

M crit  ...  (1) 


where  E*  is  the  web  material  elastic  modulus. 

The  buckling  coefficient  K depends  on  the  panel 
boundary  conditions  and  the  panel  proportions,  given  by 


♦ 


c 


c* 

a* 


(2) 


The  panels  will  be  assumed  simply  supported  at  the 

long  edges  of  the  web. 

* 

The  remaining  boundaries  are  provided  by  the  stiffeners, 
whose  influence  depends  upon  the  stiffener  bending  stiffness 
parameter 


}i 


E*I*c* 
s s 

5 — r 

E*a*  t* 


ll 


where  I * stiffener  effective  second  moment  of  area. 

8 

E_  ■ stiffener  material  elastic  modulus, 
s 

Reference  9 gives  the  buckling  coefficient  K as  a 
function  of  c and  u in  graphical  form. 

From  equation  (1),  the  skin  thickness  required  to 
prevent  buckling  due  to  the.  shear  flow  q*  is 


From  equation  (3),  stiffener  bending  stiffness  required  is 


E*  I*  - 

s 8 c* 

of,  if  stiffener  and  web  are  made  from  the  same  material 


j*  . u^a^t*3 

8 c* 


...  (5) 


By  assuming  a linear  relationship  between  stiffener 
bending  stiffness  and  stiffener  area,  a dimensionless 
stiffener  shape  factor  6 may  be  determined  for  any  given 

S 

family  of  stiffener  designs,  so  that 


I.  - ...  (6) 

Equation  (6)  would  apply,  for  instance,  to  a range  of 
thin-ifalled  stiffeners  of  constant  size,  but  with  varying 
wall  thickness. 

The  equivalent  web  thickness  attributable  to  the 
stiffeners  may  now  be  found  from  (5)  and  (6)  which  give 


C*  V*’ 


so  that  from  (4\ 


U*q* 

kO> 
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...  (7) 


An  expression  for  total  equivalent  thickness  may  now 
be  written  in  the  following  form 


e ? 


Hit] 


* R-S 


...  (8) 


t*  t*  t* 

where  te  = — = — + — 

a*  a*  a* 


...  (9). 
...  CIO) 


The  quantity  n is  the  skin  efficiency  factor  which 
will  be  defined  in  appendix  A for  the  foam  stabilised  web. 

If  the  panels  between  stiffeners  are  Unstabilised,  i.e.  without 
additional  foam,  n ■ 1. 

Using  equation  (8),  for  any  given  value  of  the  parameter 
R the  quantities  c and  y may  be  varied  to  determine  their 
optimum  values,  so  that  the  total  equivalent  thickness  is 
minimised. 

Since  the  buckling  coefficient  K is  available  as  a 
function  of  c and  p only  in  numerical  form,  the  optimisation 
process  will  inevitably  involve  numerical  interpolation. 

However,  experience  shows  equation  (8)  to  be  well-behaved  so 
that  the  required  results  may  be  obtained  numerically  with 
excellent  consistency  and  accuracy. 

A Fortran  routine  has  been  developed  which  utilises  the 
available  numerical  data  to  evaluate  equation  (8). 

For  given  values  of  load  intensity  R and  web  efficient  n 
the  programme  determines  values  of  the  equivalent  thickness 
parameter  Tg  as  a function  of  c and  y,  and  then  proceeds  to 
Identify  the  locus  of  minimum  T0  values. 

Data  generated  by  this  programme  is  shown  plotted  in 
figures  2 to  4 . A second  version  of  the  programme  is  available 
which  operates  entirely  with  dimensional  quantities.  Both 
programmes  are  described  in  Appendix  B. 

2.2  Orthotropic  plate  analysis 

The  design  procedure  described  below  is  based  on  the 
analysis  of  Timoshenko  and  Gere  , page  407,  which  assumes 
web  bending  stiffness  to  be  uniform  at  all  points  of  the  plate, 
although  having  different  values  in  direction  parallel  and 
perpendicular  to  the  plate  edges. 

Reference  11  gives  for  the  plate  critical  shear  flow,  .. 


<*?rit 


where  Dj  a 


D, 


“(DtDj*  )* 

a** 


...  (ID 


plate  bending  stiffness  perpendicular  to  stiffeners 
E*t*J 


12 


...  (12) 


plate  bending  stiffness  parallel  to  stiffeners 

...  (13) 


E* 

s s 

c* 


k 

0 

e 

b* 


buckling  coefficient  depending  on  parameters 

e.e. 

a* 

F* 


r M1 

lp2J 


(DjD?J. 

D, 


i 


...  (14) 


overall  web  dimension  perpendicular  to  stiffeners 
(■  ® in  this  case) 


J)s  ■ plate  twisting  stiffness 
G*t*> 


(15) 


When  B ■ 0,  the  following  equation  fits  the  data  given  by 

Timoshenko  in  figure  9-42  with  good  accuracy 


8 + — 
8 e 


(16) 


Equations  (12)  to  (16)  may  now  be  substituted  into  (11)  so  that, 
after  some  reorganisation,  the  following  cubic  equation  is 
found 

(TI  *)’  ♦ — (TI*  ) - if*  - 0 ...(17) 

s 16/3(l+v)  s 

I*  t 

where  I * — ; T * — U ...  (18) 

s c*t**  Q7 

Equ  ition  (17)  is  in  Cardan's  form,  x*  + rx  ♦ s = 0,  for  which 
the  following  root  exists 


X 3 


(19) 
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f.  £ 4 (s2+rli]K+  r s _ 

L 2 '4  27'  J [2  '4  27  ' J 


where 


TI 


8 


5T2 


16/3(l«-v) 
s - 12* 

8 vr 

It  may  be  demonstrated  that  for  the  range  of  values  involved  in 

the  current  problem 

.2 


r^ 

27 


s‘ 

<<  4 


so  that  (19)  simplifies  to  give 


lj  T 


(S') 


. . . (20) 


The  total  equivalent  plate  thickness  may  be  formed  by  using 
equations  (6),  (9),  (10)  and  (18)  to  give 

Te  “ T*f(IB*T)% 


which,  using  (20)  gives 


♦ £8 


...  (21) 


The  value  of  T to  give  minimum  Tg  may  be  found  from  (21) 
by  differentiation,  which  process  produces  the  following  results. 


T = p 2 /H 
so 


T - 2p2/K 

eo 


I 


so  pzR 


. . . (22) 
...  (23) 
. . . (24) 


A range  of  optimum  designs  of  equal  merit  lie  on  a 
ray  in  u-c  space  given  by 

c * p2Ry  . . . (25) 

Expressed  in  dimensional  quantities,  the  equivalent 
shear  stress  for  optimum  designs  is  given  by 


. ..  (26) 


eo 


2p*  taVj  . 


and  since  T 


so 


T . the  true  stress  in  the  web  is 
o 


2t* 

eo 


(27) 


The  remainder  of  the  y- c space  may  be  mapped  by  noting  from 
(3)  and  (18)  that 


I 


8 


so  that  from  (20) 

c * 

T - P(^)  ...  (28) 

This  may  be  substituted  into  (21)  which  may  then  be  solved  to 
give  an  equation  for  the  equivalent  thickness  contours  in  the 
following  form. 

c ‘ [(sp) 1 {(sf)  ' p,^J  • “ •••  (29) 

The  lines  are  rays  from  the  origin,  the  upper  and  lower  roots 
giving  contours  on  either  side  of  the  optimum  given  by 
equation  (25).  A typical  solution  is  illustrated  in  figure  5. 

Care  must  be  exercised  in  applying  equation  (29)  to  the 
design  process.  Orthotropic  plate  theory  has  been  shown  to  apply 
reliably  only  when  stiffeners  are  very  closely  spaced  (see 
reference  8).  Thus  equation  (29)  should  be  regarded  as  simply 
providing  a guide  to  behaviour  in  the  extreme  regions  of  the 
design  charts,  as  shown  in  figures  2 to  4. 

3.  Discussion 

Certain  consistent  features  may  be  observed  from  the  results 

based  on  equation  (8),  which  are  illustrated  in  figures  2,3,4  and  9 

/ 

It  would  appear  that  for  most  practical  purposes,  weight  reduced  S 
continuously  as  stiffener  spacing  is  reduced. 


This  process  reaches  its  limit  at  very  small  stiffener 
spacings  indicated  by  the  orthotropic  plate  solution  in 
equations  (22)  to  (25).  This  solution  appears  to  be  consistent 
with  that  based  on  the  more  accurate  analysis  which  takes 
account  of  finite  stiffener  spacing  in  the  buckling  mode. 

Numerical  results  based  on  the  more  sophisticated  analysis  are 
available  only  for  values  of  c greater  than  0.2,  and  therefore 
this  has  been  chosen  as  a point  of  demarcation  between  the 
solutions. 

The  addition  of  low  density  material  to  plate  elements 
offers  a substantial  inprovement  in  efficiency,  as  indicated 
in  figure  8.  Figure  10  indicates  that  these  improvements  may 
be  successfully  incorporated  in  shear  webs  of  the  type  under 
consideration. 

One  way  in  which  these  benefits  may  be  realised  is  to 
utlise  a higher  web  efficiency  to  allow  larger  stiffener 
spacing  without  weight  penalty,  thus  producing  a more  practical 
design. 

4 . Conclusions 

An  analysis  has  been  developed  which  leads  to  optimum 
designs  for  long  shear  webs  with  simply  supported  edges 
stabilised  v a regular  array  of  vertical  stiffners.  Stress 
levels  are  . ssumed  to  be  below  the  elastic  limit,  and  a linear 
relationship  between  stiffener  area  and  second  moment  of  area 
is  implied. 

The  effect  of  additional  low  density  material  applied  to 
the  skins  between  stiffners  has  also  been  studied. 

The  results  indicate  that  structure  weight  decreases 
continuously  as  stiffener  spacing  is  reduced  until  an  absolute 
minimum  is  achieved  at  a spacing  well  below  practical  limits. 

The  benefits  to  be  obtained  from  the  addition  of  low  density 
elements  are  found  to  be  substantial. 

Computer  routines  based  on  the  analysis  have  been  written. 
These  programmes  may  be  used  as  a basis  for  further  parametric 
studies  (SHEA)  and  fer  practical  design  purposes  (S!YEB). 


5.  Principal  notation 
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c* 
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AS 


xs 


q* 

E* 

G* 

T* 

p* 
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t 
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T 

Z 

n 

Dw 

e,e 
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Web  depth 
Stiffener  spacing 
skin  thickness 

stiffener  equivalent  skin  thickness 
total  equivalent  skin  thickness 
stiffener  area 

stiffener  second  moment  of  area 

shear  flow 

Young's  modulus 

shear  modulus 

shear  stress 

density 


«*t* 

c*/a* 

t*/a* 


c Isl 


« »q* 

a*E* 


/♦ 


t/Q 
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plate  bending  stiffness  parameter  (equation  A3) 
skin  efficiency  factor 
orthotropic  plate  bending  stiffnesses 
orthotropic  plate  parameters 
orthotropic  plate  buckling  coefficient 


skin  buckling  coefficient 

(S')" 


j 
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The  Effect  of  a Secondary  Stiffening  Layer  on  Plate  Buckling 


Consider  a thin  flat  plate  of  uniform  thickness  with 
given  boundary  conditions. 

The  edge  loading  at  which  such  a plate  will  buckle  is 
given  by  equation  (1)  which  may  be  re-written 


q* 

crit 


c* 


...  (Al) 


where  I* 


plate  second  moment  of  area  per  unit  width 
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We  may  now  consider  a similar  plate  composed  of  two 
layers  of  isotropic  material  characterised  by  suffices  1 and  2, 
the  two  materials  having  a common  value  of  Poisson's  ratio  but 
being  otherwise  different.  By  assuming  a linear  distribution 
of  bending  strain  through  the  thickness  of  the  combined  plate 
the  following  expression  may  be  derived 

ZEftf* 

(E*I*)  - — ...  (A2) 


where 


1 + 


b-i-V 


t xi 

i) 


...  (A3) 


If  equations  A2  and  A3  are  substituted  into  Al , the 
thickness  t*  of  material  1 may  ve  found  which  is  required  to 
prevent  buckling  under  a given  loading  q*  in  the  following 
form 


• • • 


(A4) 


The  weight  of  the  plate  per  unit  area  is  given  by 

**  - f PSt]*  • 


...  (A5) 


where 


i ♦ « . 21‘ 

pT  t* 


...  (A6) 


p*  * respective  material  densities 
»#  2 

lVhen  plate  2 is  absent,  t*  = 0 and  n * 1,  and  so  n,  which  is  a 
function  of 


P*  . E* 


i . fy  and  7? 

i E i 


represents  the  efficiency  of 


the  combination,  relative  to  the  single  layer  plate  of  material  1 


A t/  X 

As  — •,  n 


T -fej  - n’ 


...  (A7) 


n2  is  the  relative  efficiency  of  a plate  composed  entirely  of  material 
2.  Equation  (A7)  indicates  why  low  density  materials  are  superior 
for  stability-critical  structures,  since  for  most  metals 


E*  o p* 


.*.  from  (A7) 


l€ 


...  (A8) 


If  material  2 is  a low  density  organic  foam  and  material  1 
is  a metal,  p*  will  be  an  order  of  magnitude  less  than  p*  and 
although  usually  for  this  combination  of  material 

Ei  Pi 

1J  potential  efficiency  is  high. 


Much  of  this  gain  in  efficiency  may  be  achieved  without 
an  excessive  thickness  of  foam,  as  is  shown  in  figure  7. 


The  resulting  structure  is  simple  to  manufacture,  with 
few  detail  design  problems. 

The  above  analysis  applies  to  the  design  of  any  thin 
plate  element  which  is  stability  critical  under  any  given 
combination  of  edge  loading. 

Thus  the  results  may  be  applied  to  the  plate  elements 
between  the  stiffeners  of  the  shear  webs  which  are  the  subject 
of  this  report. 

The  substantial  benefit  to  be  obtained  from  secondary 


Appendix  B 

Computer  Programmes  for  Shear  Web  Design 

Two  For tan  programmes  are  available  SHEA  and  SWEB,  which 
are  iTased  on  equation  (8). 

SHEA 

This  programme  works  entirely  in  terms  of  non-dimensional 
quantities  and  is  suitable  either  for  terminal  operation  or 
batch  processing.  Input  requires  two  numbers  only,  the  load 
intensity  parameter  R and  skin  efficiency  factor,  n. 

Optimum  solutions  are  given  for  a range  of  stiffener 
spacings.  Since  design  continuously  improves  as  stiffener 
spacing  is  reduced,  the  designer  must  choose  from  this  list 
the  best  design  available,  which  will  be  limited  either  by  the 
smallest  acceptable  stiffener  spacing,  or  the  smallest  acceptable 
stiffener^size. 

If  required,  a complete  carpet  of  designs  may  be  output 
as  an  aid  to  parametric  studies. 

Run  times  depend  on  installation,  but  are  brief. 

• 

SWEB 

This  programme  is  similar  to  SHEA,  but  uses  dimensional 
data  for  more  direct  design  applications. 

Input  includes  shear  force,  web  depth,  Young's  Modulus, 
a sample  stiffener  area  and  corresponding  second  moment  of 
area,  and  skin  efficiency  factor. 

Output  which  is  dimensional,  is  restricted  to  optimum 
designs,  but  includes  information  on  web  stresses.  Consistent 
units  are  required  throughout. 

Listings  and/or  card  decks  are  available  from  the  author. 
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FIGURE  3 


Shear  web  equivalent  thickness  related  to  stiffener 
spacing  and  stiffness 
(ii)  R=0.5;  n«l .0. 
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FIGURE  10.  Shear  web  equivalent  thickness  related  to  skin 
efficiency.  (R  « 0.5) 
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